Statistics of photodissociation spectra: nonuniversal properties 
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We consider the two-point correlation function of the photodissociation cross section in molecules 
where the fragmentation process is indirect, passing through resonances above the dissociation 
threshold. In the limit of overlapping resonances, a formula is derived, relating this correlation 
function to the behavior of the corresponding classical system. It is shown that nonuniversal features 
of the two-point correlation function may have significant experimental manifestations. 
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The photodissociation of molecules, such as the radi- 
cals HO2 and NO2, is an indirect process comprised of 
two steps. In the first step, a photon excites the molecule 
to an energy above the dissociation threshold. Then frag- 
mentation proceeds via energy redistribution within the 
vibrational degrees of freedom, or tunneling from binding 
to unbinding energy surfaces of the adiabatic electronic 
potentialtl. A barrier, which separates quasi-stable states 
from continuum modes, hinders the immediate dissocia- 
tion of the excited molecule. On these long lived reso- 
nances, the dynamics of the system is chaotic. Therefore, 
the photodissociation cross section, exhibits a. compli- 
cated behavior as function of the photon energyH, behav- 
ior which suggests a statistical analysis of the problem. 
This appsftach has been taken, recently, by Fyodorov and 
AlhassidHa who analyzed correlations of the photodisso- 
ciatiorj-|Cross section in the framework of random matrix 
theoryH. In this framework, one is able to describe the 
universal properties of the photodissociation cross sec- 
tion, but not the individual imprints of each molecule. 
It is well known, however, that the statistical proper- 
ties of quantum chaotic systems bear a simple relation to 
the underlying classical dynamicaj. The main purpose of 
this paper is to relate the statistical characteristics of the 
photodissociation cross section to correlation functions 
associated with the classical dynamics of the molecule. 

Consider a molecule, in the ground state \g), excited by 
a light pulse to an energy above the dissociation thresh- 
old, and let H denote the Hamiltonian of the system on 
the excited electronic surface. It will be assumed that 
Ti represents an open system with several open dissocia- 
tion channels. The photodissociation cross section of the 
molecule, in the dipole approximation, is given by 



cr(e) = -?/ Ini {g\ V- 
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where e is the energy measured from the ground state 
of the molecule and e"^ — e ± iQ, V — d ■ e is the pro- 
jection of the electronic dipole moment operator of the 
molecule, d, on the polarization, e, of the absorbed light, 
and 77 = e/cheo, c being the speed of light, and cq the 
electric permitivity. 

We focus our attention on the dimensionless two-point 
correlation function: 
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where (• ■ •) denotes an energy averaging over a classically 
small energy interval centered at e which, nonetheless, in- 
cludes a large number of resonances. It will be assumed 
that the excitation energy, e, is sufficiently high such that 
the mean spacing between the vibrational modes of the 
molecule is smaller than the energy broadening due to 
the finite life time of the system in the excited states. 
This is the regime of overlapping resonances. 

Turning to the classical counterpart of our system, let 
p(x) be the initial density distribution of the excited 
molecule in phase space. Here x = (r, p) is a phase space 
point, r and p being d dimensional vectors of coordinates 
and conjugate momenta, respectively. The classical au- 
tocorrelation function of p(x) is defined as 



C(t) 



((P(x(t))p(x))) 



where x(i) is the end point of the trajectory starting at 
x(0) = X, and evolving for a time t according to the clas- 
sical equations of motion: f ~ dH/dp and p = —dH/dr. 
The averaging, ((• • ■)), is over all initial points, x, on the 
energy shell e = 'H(x) (see definition in Eq. (|l^) below). 
Another classical correlation function, closely related to 
C(t), is defined by restricting the average in (||) only to 
those points x which lie on the periodic orbits of the sys- 
tem with period r, where t > t. This function will be 
denoted by Cr{t), and its precise definition will be given 
later on (see Eq. (|l^)). 

The central result of this paper is a formula which 
expresses the two-point correlation function of the pho- 
todissociation cross section, K{u!), in terms of the clas- 
sical correlation functions, C{t), and Cr{t). Namely, 
within the semiclassical approximation, K{aj) takes the 
form 
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Here (3=1 for Hamiltonians with time reversal symme- 
try, and /3 = 2 for systems without this symmetry. It is 
assumed that H have no other discrete symmetry. 

The above results will be derived by semiclassical 
methodscl. The overlap of resonances implies that the be- 
havior of CT(e) is dominated by short time dynamics, and 
therefore the semiclassical approximation is justified. A 
semiclassical treatment of this problem will also clarify 
the role of various classical trajectories, and help inter- 
preting formulae (|^^. 

In the semiclassical limit, the Green function is approx- 
imated by a sum of two terms: G ~ Gw + Gosc- The 
first, known as the Weyl term, is a smooth function of 
the energy, e: 

{2Trhy e± - H (p, ^) ^ ' 

It is large when |r — r'| is of order of the particle wave- 
length, and therefore represents a local function in space. 
The second contribution, Gosc, is nonlocal in space but 
oscillatory in the energy. It is expressed as a sum over 
classical trajectoriesB: 

(r'|G„,,(6±)|r)=^A±e±^^-(^^'-'''-), (8) 

where Syit; r', r) is the classical action of the i^-th trajec- 
tory from r to r' with energy e, and = {A^)* is the 
corresponding amplitude. 

This decomposition of the Green function implies that 
the average, {a{e)), comes only from the Weyl contribu- 
tion, Gw- A straightforward calculation yields {(j{e)) = 
((pW)), where 

p(x)= y'dr'e^P--'(r-^|V)(V'|r + ^) (9) 

is the Wigner function of the (unnormalized) state \4>) = 
'D\g). This function is real, and upon small smearing in 
phase space yields a positive definite function which can 
be interpreted as the initial density distribution of the 
excited molecule. The microcanonical average, (( )), 
of a general function /(x) is defined as 

((/(x))) = J S{e 7i(x))/(x). (10) 

Notice that / and ((/)) do not have the same dimensions. 
To simplify our notations, from now on we work in units 
where h = 1. 

The correlator of cr(e) can be written in the form: 




xRe(G(e--|;r2,ri)G(e+ + |;r4,r3)), (11) 

where p(r, r') = (r|?/'} (■0|r'). Considering the connected 
part of the correlator, only the oscillatory terms of 



the Green functions, Gosc, contribute. Since these are 
strongly oscillating functions of the coordinates, the main 
contribution to the integral (|l^) comes from two types of 
orbits: (a) general orbits with initial and final points far 
apart, but r2 ~ rs and ri ~ r4. (b) Returning orbits 
with ri ^ r2 and rs ~ r^. The corresponding contribu- 
tions to (|l|) will be denoted by la and respectively. 
When Ti. is also time reversal symmetric, additional con- 
tribution, equal to la, comes from orbits with ri ^ rs 
and r2 ~ r4. 

The calculation of la involves two main approxima- 
tions. One is to expand the actions of the Green functions 
(^ to the first order in — T2, and — ri. For this 
purpose we use the relation: Si,{e±uj/2; r' + 5r' ,r + 6r) ~ 
S,y{e; r', r)±t^uj/2 + p'j^ ■ Sr' — p^, ■ Sr, where t^, is the time 
of the orbit, Pi, is the initial momentum of the trajectory 
(at point r), and p'^^ is the final momentum (at point r'). 
The second approximation is to neglect quantum inter- 
ference effects (weak localization corrections) and employ 
the "diagonal approximation" . In this approximation, 
the double sum on trajectories in the product of Green 
functions, Goscie+)Goscie-) ^ Y.^.^. A,Al,e'S--'S., ^ 
replaced by a single sum in which one keeps only pairs of 
trajectories, {v,v'), related by symmetry. Noticing that 
/3(ri,r2)p(r3,r4) = p(ri, r4)p(r3, r2), and integrating on 
the coordinate differences — V2, and — ri, one thus 
obtains: 

la ^ ^Re j dvdv' ^ V(xL)p(x.)e'*^-, (12) 

where = (r,p,y) and x^ = (r',p'j,) are the initial and 
final phase-space coordinates of the i/-th trajectory. |_. 
Now, to evaluate the sum over z/, we use the sum rulell3: 

(2^)'^-i5]|A,p5(x:,,x,,t,)= (13) 
/ dtdp'dp g(x', X, t)5{e - H(x))5(x' - x(t)), 

JQ 

where g{x',ji,t) is a general function of the time t, the 
initial x, and the final x' phase space points. x(t) is the 
phase space coordinate, after time t, of a trajectory start- 
ing from x(0) = X. Eq. (|l^) can be proved by performing 
the integrals on its right hand side. For this purpose it 
is convenient to use a local coordinate system with one 
cooridinate parallel to the orbit, and d — I perpendicular 
to itt3. Substituting (|l|) in (|l2|), and integrating over X 
yields 

la ^ ji^^^J/i Jd^ e**-<5(e-W(x))p(x(t))p(x). 

(14) 

Using the definition of microcanonical average (M) and 
the correlation function (^, one obtains KiJlu) = 



2 



2/a/((p))^/3, where (5 is the symmetry factor associated 
with the additional contribution in the case of time re- 
versal symmetry. 

To evaluate the term associated with returning tra- 
jectories, we expand the actions of the Green functions 
in to first order in V2 — ri and r4 — r3. Integrat- 

ing over these coordinate differences, the result takes the 
form h = ifRe{I{Lo)r{Lo)} /27r2, where 



dv y ^ A„p(xt.)e 



iS„(e+^;r,r) 



(15) 



Here x^ = (r, p^,) is a phase space point correspond- 
ing to the v-th. returning trajectory at point r, where 
Pi/ — {v'u + Pi')/2 is the average of the initial, Pi,, and 
final, p'j,, momentum of the trajectory. The integral over 
r is now performed in the stationary phase approxima- 
tion. The stationary phase condition, dSv{e;r,r)/dr = 
P'u = implies that the initial and final momentum 
of the orbit are equal, therefore, the main contribution 
to th&iptegral ( |T5| ) comes from the vicinities of periodic 
orbitsB. The result can be expressed as a sum over peri- 
odic orbits: I{lu) = J^p^Ppe'^"^'^^'*"'^^'^, where Ap, tp 
and Sp are the amplitude, period, and action of the p-th 
periodic orbit respectively, pp — § dtp{-x.p{t)) is the time 
integral of p(x) along the orbit Xp(t). 

Next, we analyze the product I{lu)I*{lu), which in the 
diagonal approximation is. 
Following Eckhardt et al.ll 
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\t dt'p{^p{t - U)p{-Kp{t + U) = 



dt'{p{^{t'))p{^))^^ 



where ((• • ■))p — ip ^ / dt{- ■ ■) denotes the time aver- 
age along the p-th. periodic orbit. Thus {I{lu)I*{lu)) — 
{2/f3)J2tp\Ap\'^e^*p'^{pix{t'))p{x))p. The sum over the 
periodic orbits can be converted into an integral over 
time using the relation 



(16) 



t<tp<t+St 



= St JdxS[e^n{x)]5[x\\ 



where X|| is a phase space coordinate on the energy shell 
7i(x) ~ e, while X|| {t) is the position of the system af- 
ter time t, starting from X||(0) — X||. Using (|l6|), the 
sum over periodic orbits is converted into a time integral 
which yields formula (|^) for K2{uj), with the correlation 
function: 



Crit) = jd^ 



P(x(0)p(x) 



<5[e-W(x)]<5[x,|-x,|(r)]. (17) 



The identification of the classical trajectories which 
underly the main contribution to the correlator K(lu), 
clarifies the origin of fluctuations in cr(e). These have 
two sources: the wave functions of the system, and the 
density of states. Ki{uj), associated with general trajec- 
tories, is related to fluctuations in the wave functions, 
while K2{i-j), coming from the periodic orbits of the sys- 
tem, is related to fluctuations in the density of states. 

To proceed further, one need to characterize the be- 
havior of the correlation functions C(t) and Cr{t). For 
this purpose we assume that the classical dynamics of 
the system is characterized by two time scales which are 
well separated. Such a situation appears, for example, 
in almost close chaotic systems. There the long time 
scale is the typical time for dissociation, while the short 
time scale is the time which takes for a classical den- 
sity distribution to relax to the ergodic distribution on 
the energy shell. We denote by Ho the closed Hamilto- 
nian which controls the dynamics of the system for time 
shorter than the dissociation time. Under these assump- 
tions C{t) ~ ACtir) ~ e^'''l*l/A, where 7 is the dissocia- 
tion rate, and A is the mean spacing between resonances, 
i.e. (27r;i)''/A = J dxS{e - ■Wo(x)). We assume A to be 
approximately constant within the interval of energy av- 
eraging, and from now on work in units where A = 1. If 
A is not constant, a{e) should be unfolded appropriately. 

An infinite separation between the time scales men- 
tionedpjibove corresponds to the limit of random matrix 
theoryu. In this case K{uj) reduces to 
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(3tt \u!^ + 72 27r (72 + uj^y 



(18) 



The first and the second terms of this formula come from 
Ki{uj) and K2{uj) respectively. As can be seen, these 
components correspond to the leading and the sublead- 
ing terms of an expansion in I/7. 

Formula (|l^), describing the correlations in the univeiu 
sal regime, was first derived bj« Fyodorov and Alhassidu 
using the nonlinear cr-modelt3. Here we confirm their 
conjecture that for overlapping resonances (7 > 1) K{llj) 
can be derived by semiclassical methods. 

Yet, the range of validity of formulae (^||) goes far be- 
yond the universal regime. They account also for system 
specific contributions to K{uj), which might be of the 
same order of the universal result (|l^). To be concrete, 
we focus on the nonuniversal corrections due to the lead- 
ing term, Ki{u!), and consider a representative behavior 
of decay of correlations in open chaotic systems: 



C(t) ~ e'^l*! -hQ^e-^=l*lcos(w2^)• 



(19) 



Here 72 is the real part of the second Ruelle resonanc 
of the classical system, while 1^2 is its imaginary part, a 
is a constant of order unity which depends on details of 
the system and the initial state (^) . The situation where 
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'^2 7^ appears usually when a specific short periodic or- 
bit has strong influence on the dynamics of the system, 
since every typical trajectory stays in its vicinity for a 
long time. Thus, 0^2 ~ iir/t* , where t* is the period of 
the orbit. Substitution of @ in (|-|) yields: 
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FIG. 1. The two-point correlation function of the photodis- 
sociation cross section for systems with time reversal symme- 
try (/3 = 1). The solid line is the universal result Ku['y,oj) 
(Eq. 18) for the case 7 = 1. The dashed and the dotted lines 
correspond to the correlation function, K{lo), which includes 
the nonuniversal contribution, in the approximation of Eq. 
(20). In both cases 7 = 1, 72 = 3 and q = 1, but, the dashed 
line corresponds to ulJ2 = 0, while the dotted line to U2 ~ 6. 

This formula is plotted for various parameters in Fig. 1. 
The solid line represents the universal form (p^. The 
dashed line, corresponding to UJ2 = 0, shows the typical 
behavior of systems where the decay of correlations is of 
diffusive nature. In this case the main deviation from 
universality appears as an increase of correlations near 
the peak of K{lu), at w = 0. The dotted line corresponds 
to a nonzero value of U!2, and characterizes the behavior 
of ballistic systems. Here the nonuniversal contribution 
is located in the tail of K{uj), near uj — uj2, where the 
universal term ( |l8| ) is already negligible. These plots 
demonstrate the significance of system specific contribu- 
tions to K{u!). 

In conclusion, we derived the semiclassical relation be- 
tween the two point correlation function, K{uj), and the 
classical correlation functions C{t), and Cr{t). The first 
quantity characterizes the quantum mechanical process 



of photodissociation, while G{t), and Cr(<) characterizes 
the dynamics of the classical counterpart of the system. 
In contrast with the two-point correlation function of the 
density of states, where periodic orbits play the major 
role, here we have found that the main contribution to 
K{ljj) comes from general orbits. Moreover, the nonuni- 
versal contributions to the correlator, K{llj), can be of the 
same order of the universal result of random matrix the- 
ory (the relative magnitude being proportional to 7/72). 
These results can be used in order to analyze the exper- 
imental data of complex molecules, extract information 
regarding their classical dynamics, and construct effec- 
tive models for these molecules. 
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